Universal non stationary dynamics at the depinning transition 
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We study the non-stationary dynamics of an elastic interface in a disordered medium at the depin- 
ning transition. We compute the two-time response and correlation functions, found to be universal 
and characterized by two independent critical exponents. We find a good agreement between two- 
loop Functional Renormalization Group calculations and molecular dynamics simulations for the 
scaling forms, and for the response aging exponent On. We also describe a dynamical dimensional 
crossover, observed at long times in the relaxation of a finite system. Our results are relevant for 
the non-steady driven dynamics of domain walls in ferromagnetic films and contact lines in wetting. 
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The universal glassy properties that emerge from the 
frustrating competition between elasticity and disorder 
are relevant for many experimental systems, such as in- 
terfaces describing magnetic [H, [3, 9 and ferroelec- 
tric Q domain walls, contact lines of fluids 0, [1] and 
fracture Disorder leads to pinning, affecting in 

a dramatic way their dynamical properties. In partic- 
ular, when driven by an external force / at zero tem- 
perature, disorder leads to a depinning transition at a 
threshold value / = /c, below which the interface is im- 
mobile, and above which steady-state motion sets in. For 
/ ^ /c it has been fruitful to regard the depinning tran- 
sition as a critical phenomenon, with the mean veloc- 
ity V as an order parameter, v ^ {f — fc)^ , and with 
a characteristic length ^ playing the role of the diver- 
gent correlation length ^ ~ (/ — fc)'", P and v being 
universal exponents Near the critical point, how- 

ever, the time needed to reach such a non- equilibrium 
steady-state can be very long, since the memory of the 
initial condition persists for length scales larger than a 
growing correlation length i(t) ^ t^^^, with z the dy- 
namical exponent 1^ Il3l |. Being only limited by the 



divergent steady correlation length ^ or the system size 
L, the resulting non-steady critical regime is macroscopi- 
cally large, t ^ £^^,L^. It is thus relevant for experimental 
protocols. Analogously to non-driven systems relaxing 
to their critical equilibrium states IJ, , we show here 
that the transient dynamics of a driven disordered system 
displays interesting, though different, universal features. 



Dynamical properties are characterized by two-time 
t,tw response and correlation functions, describing the 
time evolution of the system as a function of its "age" 
or waiting-time tw from a given initial condition at 
t = 0. Progress was achieved in the understanding of 
the steady-state, where these functions depend only on 
t — tw Functional Renormalization Group (FRG) cal- 
culations 3, 12, 18, 3 allowed to compute the critical 
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FIG. 1: Scaling of the local Ptt^" and the center of mass Ptt^" 
integrated linear response functions of an elastic string re- 
laxing at its depinning threshold (they have been shifted for 
clarity). The new critical exponent Or = —0.6 ± 0.05 numer- 
ically obtained is in good agreement with our two-loop FRG 
prediction Or ~ —0.64 (|12p . Here z — 1.5 is the dynamical 
exponent, estimated numerically. Inset: non-scaled response 
data for the center of mass. 



exponents describing different universality classes, and 
powerful algorithms were developed to elucidate the low 
temperature dynamical phase diagram f20j. In contrast, 
little is known about the more difficult transient regime, 
where the time-translation invariance is broken. Yet the 
first steps in that direction have unveiled rich and uni- 
versal behaviors, including slow dynamics [13J and aging 
properties characterized by new exponents [12| . 

Here we focus on elastic interfaces of dimension d {d = 
1 for an elastic line) parameterized by a scalar field Ux,t 
describing their position in a d-\~ 1-dimensionaI disordered 
medium. The driven ovcrdampcd dynamics of this model 
system obeys the equation of motion 
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where rj is the friction coefficient, c the elastic constant, 
and F{x,u) a quenched random pinning force with dis- 
ordered averaged correlations F{x,u)F(x' ,u') — A(u — 
u')S'^{x — x'). Under an apphed force /, the velocity is 
V = L^'^ j d'^xdtUx^t- In this paper we consider a flat 
initial configuration Ux^t=Q — but our results hold for 
any short ranged correlated initial conditions. Denoting 
iiq^t the spatial Fourier transform of Ux,t, we focus on the 
linear response Ti-tt^ to a small external field h-qt„ and 
the correlation function C?, : 



The main result of this Letter is to establish, both via 
numerical calculations and additional analytical work, 
that these central observables take the scaling form: 

= {t/t^f''q'-''FR[q'-{t-t^),t/t^], (3) 
Ctt^ = q-^'+^^\t/t^Y^-^Fc[q^{t-t^),t/t^]A^) 

with Or and 9c two new universal critical exponents, i.e. 
independent of the usual depinning exponents. These 
are defined such that FR_c(yi,?/2 ^ oo) fR,c{yi) 
for fixed yi. In Eq. ([3]), dH), Fr^q are universal scal- 
ing functions (up to a non universal amplitude) and ( 
the roughness exponent. In the limit yi — > 0, 2/2 fixed, 
one finds ^^(2/1,2/2) ~ 2/f "^■'^''gi?(y2), -^0(2/1,2/2) ~ 
^(d+2Q/z ^^^^^^^ j ^ defined q — ^ limit. These 
scaling forms were predicted in Ref. |I3] based on a 
one-loop FRG calculation. One may question however 
whether this lowest order in the d — A — e dimensional 
expansion is accurate enough to describe interfaces of ex- 
perimental interest d — 1,2. In addition, no prediction 
for 9c was obtained. Here we firmly establish that the 
above scaling forms hold and we provide a reliable de- 
termination of Or and Oc vn d — 1. We also perform a 
two-loop FRG calculation, as is known to be required for 
a consistent theory of depinning [19]. 

Most of the numerical studies of the transient dynamics 
have focused so far on one-time quantities which can be 
obtained from C^^, and nl^, in Eq. ([pand JH). The 
structure factor Sq{t) = Cl^ was found [l3| to behave as: 

Sq{t) ^ CI, ^ g-('^+2C)F[cz£(t)] , l{t) ~ t^l^- , (5) 

where F(y) ~ &^ , a constant, for y ^ 1 and F(y) ^ 
yd+2C, fpj. y The relaxational dynamics is thus dic- 

tated by a single growing length, separating the small, 
steady-state equilibrated scales, from the large ones re- 
taining a long-time memory of the initial condition. Eq. 
([5]) is obtained from (|4]) in the limit t {i.e. 2/2 1, 

2/1 — *■ 0) with q^'t {i.e. 2/12/2/(2/2 — 1)) fixed. The analogy 
with standard critical phenomena suggests for the ve- 
locity, the scaling form, v{tj) ^ h-P/''G[h-H,b^/''{f - 
fc),b~^L] where b is an arbitrary rescaling factor, nu- 
merically verified in Ref. [l3|. For f = fc and t <^ it 



implies that v{t) cc t ^/^^ and also that: 

dv{t)/df oc At^^-'^/' , (6) 

where we used the exact relations /3 — v{z — Q) and v — 
1/(2 - C), from statistical tilt symmetry (STS) [13]. We 
now check that this scaling of one time observables for 
/ = fc is consistent with the two time scaling Eq. (HI. 
Indeed Eq. ([6]) results by combining the exact relation [l8| 
d.y{t) ^ ds dtTZf^^ and the hmit g ^ of Eq. §1) 

7^«=° ^ (t - {t/t^f- gR{t/t^) , (7) 

with gR{x) oc c^* a constant for a; 3> 1. 

Let us now focus on two time quantities. To check 
Eqs. ([3]) and (|4]) we have performed numerical simula- 
tions of Eq. Il]) in the case of elastic lines, d = 1, experi- 
mentally relevant for many two dimensional systems, e.g. 
films. To study the non-stationary dynamics at f = fc 
we discretize Eq. ([T]) in the x direction, Ux,t = Ui{t), with 
i ~ 0, L — 1, and use the method described in Ref. 
We start a.t t — with a flat configuration, Ui{t = 0) = 0, 
and monitor correlation and resp onse functions at the ex- 
act sample critical force fc [2l|. Numerically, it is more 
convenient to work with the local integrated response 
= /q™ IqT^ts (where denotes the integral 
over the first Brillouin zone) and zero-mode integrated 
response " = /p " dsnl'. From Eq. © we predict, 

(i-t.)l/^ WJ' {t^t^fl^ WJ ' ^ ' 

where both h{y) and h{y) behave as y^^+^i^ for y 00. 
To implement the local (zero mode) response we define 
the observable 'Wi{t) — Ui{t) — Ucm{t) {'Wi{t) ~ Ui{t)) 
where Ucm{t) = 7; X^i ^i(^) and then compute, 

p^=°(^=°) = lim }-Y.[n,.,{t)-wf{t)]a,a-\ (9) 

i 

where wf{t) is the solution of Eq. ([T]) with uf (0) — 
Mi(0) = and an additional perturbative force aai9{tw — 
t). We take random numbers (7^ = ±1 uncorrelated from 
site to site for computing Ptt^, and CFi — (Tq for com- 
puting Ptt"". The value of a is chosen small enough to 
guarantee linear response . In Fig. [1] we show the 
numerical results for pf[^ and ptt^ ^ for L = 2048 aver- 
aged over 10000 disorder realizations. We see that the 
predicted scaling forms, Eq. (I8|), describe well the data. 
For t/t^, ^ 1 we observe a well developed power law be- 
havior with an aging exponent 9r = —0.6 ±0.05 which is 
indistinguishable for both responses, Ptt~f{t — tyj)~^/^ ~ 
p^=°(i - t„)"^/'' (i/*u.)"^+^", as predicted. How does 
this numerical estimate for 9r compare with the previ- 
ous FRG approach of Ref. [12]? The one loop result 
for 9r = + 0(e2), setting e = 3 gives 9r « -1/3. 
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Incidentally, up to one loop order, one finds the rela- 
tion 9ji = (z — 2)/z, which, using the numerical estimate 
z = 1.5 [l^, yields again On « —1/3. Although it goes 
in the right direction, it is still far from our numerical 
result. To see whether the FRG predictions can be im- 
proved we have computed 'R-lt'S *° ^°°P order 24 1 . 
The starting point is Eq. ([1]). Response and correlations 
are then obtained from the standard dynamical (disorder 
averaged) Martin-Siggia-Rose action S which reads here 



S = I iUqt[(q + dt)Stt' + ^tt']U-qt' 

't>t'>0 Jq 



iUxtiUxt' ^{Uxt - Uxt') 



(10) 



x,t>t'>0 



where A(m) is the force- force correlator and is the 
self-energy. As a result of the covariance of the action 
under STS [l3| the self-energy has the structure T.tt' — 
'^tt' — Stt' /o dtiT,tti ■ It was computed to one loop in Ref. 
[l3 | and at two loop order the perturbation theory leads 
to diagrams similar to the one contributing the dynamical 
exponent z, as depicted in Fig. 10 of Ref. [l9|, with 
the constraint that here, the time variables are positive. 



The (bare) response function TZ^^, = {iUqtU-qt' 
computed from the exact identity 



)s is then 



t'<ti<t2<t 



f>q yi jpg 



Jt'<ti<t J0<t2<ti 



(11) 



where R^^, — 6{t — t')e * is the response in the 
absence of disorder. Reexpressing in terms of A(u), cor- 
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FIG. 2: Scaling of the local autocorrelation function. The 
aging exponent is 9c ~ —1.5 ± 0.05, different from 6r. Inset: 
non-rescaled data. 



rected to the same order, and using the FRG fixed point 
equation, one explicitly shows that it has the scaling form 
as in Eq. ([7]). We then find that no new independent di- 
vergence occurs in t/t' at this order, hence that to two 



loop accuracy the relation 9r = {z-'2)/z + 0{e^) contin- 
ues to hold. Our numerical result however indicates that 
this relation cannot hold to all orders in e, i.e one must 
have On ^ {z — 2)/z, implying that 0^ is indeed a new 
independent exponent. One way to understand the FRG 
result is then to rewrite more explicitly. 
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108 648/ ^ ' 



1627\/2 

0.1111...e- 0.03395. ..e2 + 0(e3) 



(12) 



with 7 = 0.54822... We note that if we set e = 3 in 
that expression (fT2|) assuming the 0(e'^) to be small, we 
obtain On — —0.64..., very close to the the numerical 
value. Hence we conclude that although corrections to 3 
loop and higher to 0fj — {z—2)/z must be large, in 9b. they 
must be small. This provides one way of interpreting our 
results, and motivation for future analytical work. 

We have also checked numerically the scaling form for 
the correlation function in Eq. (jj]). It is more conve- 
nient to compute the autocorrelation function Cft^° = 
Ux,tUx,t^ obtained from C^^^ by integration over Fourier 
modes Cf^° = J^C^f^. From Eq. (ji]), one expects 

Cft=o - (i - ty,f'^/''h{t/ty,) where h{y) oc 
for large y. To check this, we compute numerically 
^■Wi{t)wi{t') . In Fig. [5] we show a plot 



x=0 



c, 

of {t - t^)-2C/^c^^=o foj. 2, ^ 2048 obtained by averaging 
over 10000 samples, which is in good agreement with the 
predicted scaling. For t/t^, » 1, we see a power- law be- 
havior ■v^ith a non-equilibrium exponent dp = — 1.5±0.05. 
At variance with pure critical dynamics [l5| , one obtains 
that Be and Or are different exponents. Such behavior 
was observed in other disordered elastic systems (though 
relaxing at equilibrium) [23|. The determination of 9c 
via the FRG requires a computation to order O(e^) of 
Cli^ and remains a challenge. 

So far we have analyzed the situation where i{t) < L. 
Finite-size effects can be observed in experiments how- 
ever, as shown recently for domain walls in magnetic 
nanowires NL and vortex lattices in micron-sized super- 
conductors [25[. The finite-size crossover manifests in 
the velocity as v{t)L>^/'' = f{t/L'') and for t > L"" the in- 
terface losts memory of the initial condition [l3|. What 
happens next? Remarkably, a detailed analysis of the 
dynamics for t > reveals that different initial condi- 
tions in the same sample evolve collapsing into a unique, 
sample-dependent, reparametrized velocity v{ucm) (as a 
function of the center of mass position), before stopping 
at the critical configuration, as shown in Fig.[3I^a-c). This 
behaviour can be seen as a direct consequence of the 
Middleton theorems, which assures the convergence to 
a unique periodic attractor for f > fc [1^. For t ^ 
we can thus describe the velocity by an effective equation 
of motion for a particle, iit — v{ut, f) with ut = Ucm{t). 

Before exploring further the consequences, let us note 
that the d = problem of a particle is a solvable 
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FIG. 3: (a-c) Three difFerent initial conditions evolving at 
fc in the same sample coalesce into a unique reparametrized 
velocity function v{ucm{t)) after a typical time t ^ L'' , before 
sloping at the critical configuration, Uc = Wcm(oo). (d) The 
steady-state structure of v{ucm) around Uc for different forces 
5f = f — fc ^ 0^ is a well defined parabola with a positive 
(negative) shift proportional to 5f in the position (velocity) 
axis. Inset: unshifted data. One finds for L = 32 and 5f/fc ~ 
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limit, interesting per se, as the (sample dependent) re- 
sponse function is given exactly by Rtt^ = Oit — 
tw)v{ut, f)/v{ut^, f). In the usual model r]v{ut, f) = 
F{ut) + /, near the threshold, 5f = f — fc ^ fc, the 
particle spends most time near the zero force point (set 
to be at u = 0). For a smooth force field we can write 
Sf which yields v{t) ^ t~'^ and 
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Sf>0 
Sf<0 



(13) 



for (5/ — + 0. Hence 9ji{d = 0) = —2, consistent with a 
monotonic dependence of 9fj with d, and f3{d = 0) = 1/2. 

Next, we confirm that the results of the d = model 
are relevant for the interface for £{t) > L. As shown in 
Fig. [Sji, we checked that for interfaces of sizes L = 32 
and small Sf fc the reparametrized velocity has a 
nice parabolic shape near the zero force point, v{u, f) = 
7(u -I- aSfY -I- 6(5/ where the constants b,j> (their 
size dependence will be studied elsewhere '24j). a being 
found irrelevant, this result is consistent with the steady- 
state value (3—1/2 found for the interface in the regime 
Sf^"^ » L [2^, and predicts a crossover to an effective 
= —2 in the fixed L, large t limit for the interface. 

For modeling contact lines 0] the term V'^Ux,t 
in Eq.ilT]) is replaced by a long range elastic force 
J^,[u{x' , t) — u{x, t)]/\x—x'\'^ . In this case, using the same 
numerical method, we confirm the scaling forms ([3]) (re- 



placing q^~'^ q^~^) and measure Be = —1.2 ± 0.1 and 
On = —0.5 ±0.1. The same analysis leading to (fT2| gives 
9ji — —0.22 to one loop and Or = —0.38 to two loop. 

To conclude we have confirmed numerically the scal- 
ing forms for non stationary dynamics at depinning, for 
model systems of experimental relevance. The exponent 
6ii is found in reasonable agreement with FRG predic- 
tions. An interesting dimensional crossover was found at 
large time. We hope this motivates new experiments, e.g 
in magnets and wetting. 
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